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ABSTRACT 

This  paper  examines  (linear)  secondary  instabilities  in  compressible  boundary  layers  at 
Mach  numbers  M ^  =  0,  0.8,  1.6  and  4.5.  We  find  that  there  is  a  broad-band  of  highly 
unstable  3-d  secondary  disturbances  whose  growth  rates  increase  with  increasing  primary 
wave  amplitude.  At  Mx  <  1.6,  fundamental  resonance  dominates  at  relatively  high  (2-d) 
primary  disturbance  amplitude,  while  subharmonic  resonance  is  characterized  by  a  low  (2-d) 
primary  amplitude.  At  =  4.5,  the  subharmonic  instability  which  arises  from  the  second 
mode  disturbance  is  the  strongest  type  of  secondary  instability. 

The  influence  of  the  inclination,  6 ,  of  the  primary  wave  with  respect  to  the  mean  flow 
direction  on  secondary  instability  is  investigated  at  =  1.6  for  small  to  moderate  values 
of  6.  It  is  found  that  the  strongest  fundamental  instability  occurs  when  the  primary  wave 
is  inclined  at  10°  to  the  mean  flow  direction,  although  a  2-d  primary  mode  yields  the  most 
amplified  subharmonic.  The  subharmonic  instability  at  a  high  value  of  6  (namely,  6  =  45°) 
is  also  discussed. 

Finally,  a  subset  of  the  secondary  instability  results  are  compared  against  direct  numerical 
simulations. 


*  Research  supported  by  the  National  Aeronautics  and  Space  Administration  under  con¬ 
tract  No.  NAS1-18605  while  resident  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


I.  Introduction 


In  view  of  the  renewed  interest  in  high  speed  technology  such  as  the  National  Aero- 
Space  Plane  project,  a  good  understanding  of  boundary-layer  transition  in  supersonic  and 
hypersonic  flows  is  crucial  to  the  design  of  efficient  aerodynamic  vehicles.  For  example,  an 
accurate  estimate  of  aerodynamic  drag  and  surface  heating  for  supersonic  aircraft  requires 
a  good  knowledge  of  transition  physics. 

Although  little  is  known,  either  experimentally  or  theoretically,  on  the  paths  to  tran¬ 
sition  in  compressible  boundary  layers,  it  is  well  known  that  the  early  stages  of  transition 
in  incompressible  boundary  layers  often  occur  via  (i)  primary  and  (ii)  secondary  instabil¬ 
ity  mechanisms.  Typically,  when  the  amplitude  of  the  primary  wave  (or  interchangeably, 
TS  wave)  becomes  sufficiently  large  (as  it  propagates  downstream),  the  flow  becomes  very 
unstable  to  three-dimensional  (3-d)  secondary  disturbances.  These  3-d  disturbances  are  re¬ 
sponsible  for  three  known  laminar-turbulent  transition  routes.  The  first  is  the  fundamental 
(K)  type  breakdown  in  which  the  lambda-shaped  vortices  are  aligned  along  their  peaks  in 
the  streamwise  direction,  repeating  every  TS  wavelength.  The  second  is  called  the  subhar¬ 
monic  (H)  type  breakdown:  the  lambda  vortices  are  staggered,  repeating  at  a  streamwise 
distance  equal  to  twice  the  TS  wavelength.  The  third  type  of  breakdown  is  called  detuned: 
the  spacing  of  the  lambda  vortices  lies  between  that  of  the  K  type  and  the  H  type.  The  K 
type  breakdown  often  occurs  at  a  higher  TS  amplitude  than  does  the  H-type.  Herbert  et.  al. 
[7]  have  shown  that  the  paths  to  these  breakdowns  can  be  modeled  by  a  secondary  stability 
theory  (which  studies  the  stability  of  a  spatially  periodic  base  flow);  the  K,  H,  and  detuned 
types  are  initiated  respectively  from  the  3-d  fundamental  modes,  the  3-d  subharmonic  modes, 
and  the  combination  modes. 

The  question  we  address  is  whether  the  path  to  transition  via  the  sequence  (i)  primary 
and  (ii)  secondary  instability  mechanisms  for  compressible  flows  is  viable.  We  attempt  to 
answer  this  question,  at  least  partially,  by  a  parametric  study  of  compressible  secondary 
instability  mechanisms. 

Recently,  El-hady  (1989)  [6]  and  Nayfeh  (1989)  [11]  investigated  the  subharmonic  sec¬ 
ondary  instability  of  compressible  boundary  layers  for  spatially  developing  disturbances. 
While  El-hady  studied  the  effect  of  compressibility  from  M <*,  =  0  to  M ^  =  2.2,  Nayfeh  ex¬ 
tended  the  Mach  number  range  to  4.5.  El-hady  found  that  the  local  effect  of  compressibility 
at  a  fixed  Reynolds  number  and  frequency  is  either  stabilizing  or  destabilizing,  depending 
on  the  amplitude  of  the  primary  wave  and  on  the  spanwise  wavenumber  of  the  secondary 
disturbance.  However,  when  the  primary  growth  and  the  diffusive  growth  of  the  mean  flow 
are  taken  into  account,  he  found  that  compressibility  has  an  overall  stabilizing  influence  on 
the  subharmonic  instability  modes.  Nayfeh  found  that  increasing  Mach  number  stabilizes 
the  most  unstable  subharmonic  wave.  It  should  be  pointed  out  that  their  analysis  is  based 
on  a  2-d  primary  wave,  although  for  compressible  flows,  the  most  unstable  primary  wave  can 
be  2-d  or  oblique  ( i . e. ,  the  propagation  direction  of  the  primary  wave  phase  fronts  is  inclined 
to  the  mean  flow  direction).  In  fact,  for  supersonic  Mach  numbers  less  than  3.0,  the  most 
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unstable  mode  is  always  oblique.  Clearly,  for  these  Mach  numbers,  the  effect  of  the  primary 
wavevector  angle  on  secondary  instability  mechanisms  requires  further  study.  In  this  paper, 
results  are  presented  for  both  2-d  and  oblique  primary  waves. 

In  addition  to  providing  physical  insight  into  the  physics  of  the  transition  process,  sec¬ 
ondary  instability  theory  has  several  potential  practical  objectives.  The  first  is  to  refine 
the  so-called  eN  method,  which  is  the  most  common  and  useful  engineering  tool  used  for 
transition  prediction.  The  current  e'v  method  is  based  on  the  primary  disturbance  which 
has  the  largest  total  streamwise  growth  from  the  lower  branch  neutral  point  to  the  location 
for  the  onset  of  transition.  Bushnell  [1]  pointed  out  that  a  more  refined  eN  method  should 
incorporate  disturbances  from  primary,  secondary  and  possibly  othci  highci  instability  the¬ 
ories. 

The  eN  method  has  been  quite  successful  for  transition  prediction  in  incompressible 
flows  because  of  the  following  two  reasons.  First,  the  region  from  the  incipient  growth  of  TS 
disturbance  to  the  onset  of  transition  is  dominated  by  the  TS  disturbance.  Secondly,  the 
primary  disturbance  which  has  the  largest  integrated  growth  leads  to  an  explosive  growth  of 
3-d  disturbances  which  very  rapidly  yields  to  transition.  However,  whether  the  current  eN 
method  will  work  in  compressible  flows  remain  uncertain.  Secondary  instability  theory  can 
shed  light  on  this  uncertainty. 

Secondary  instability  theory  also  provides  a  much  cheaper  avenue  for  exploratory  studies 
than  do  experiments  or  direct  numerical  simulations  (DNS).  Further,  it  provides  a  means 
to  validate  DNS  (and  vice  versa).  Relevant  parameters  (e.g.  wavenumbers,  frequencies, 
Re,  etc)  and  the  disturbance  structure  obtained  from  the  theory  can  be  used  as  the  initial 
conditions  for  DNS,  which  can  substantially  reduce  the  computer  time  required  to  simulate 
strongly  nonlinear  phenomena. 

Finally,  if  compressible  secondary  instability  can  indeed  predict  the  correct  space  and 
time  scales  and  the  structures  of  the  flow  in  the  early  stages  of  transition,  then  it  can  shed 
some  light  on  the  modeling  of  compressible  transition.  Since  no  compressible  data  (from 
experiments  or  DNS)  for  transitional  Reynolds  stresses  is  available,  this  work  may  provide 
a  blueprint  for  future  studies  and  experiments. 

The  paper  is  organized  as  follows.  Section  2  contains  the  governing  equations.  The 
primary  instability  for  compressible  boundary  layers  is  briefly  reviewed  in  section  3.  The 
secondary  instability  theory  and  its  simplifying  assumptions  are  formulated  in  section  4,  and 
in  section  5,  the  numerical  methods  are  outlined.  Section  6  contains  the  secondary  instability 
results  and  the  comparisons  with  direct  numerical  simulations.  Finally,  some  conclusions  are 
drawn  in  section  7. 
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II.  Governing  Equations 


A  schematic  of  the  boundary  layer  flows  is  sketched  in  Fig.  1.  In  the  chosen  coordinate 
system,  x,  y  and  z  denote  the  distance  in  the  streamwise,  spanwise  and  normal  directions, 
respectively.  We  denote  density,  temperature  and  pressure  respectively  by  p,  T  and  p,  and 
the  velocity  vector  by  u.  Lengths  are  non-dimensionalized  with  respect  to  the  displacement 
thickness  5*.  The  velocity,  temperature,  viscosity,  and  density  are  normalized  with  respect 
to  their  free-stream  values  (indicated  by  a  subscript  oo).  Pressure  is  normalized  with  respect 
to  pIqU^.  The  equations  governing  the  flow  of  a  viscous  compressible  ideal  gas  are  the 
compressible  Navier-Stokes  equations  and  the  thermodynamic  equation  of  state.  With  the 
aforementioned  normalizations,  these  equations  in  dimensionless  form  are 


^  +  V.(^u)  =  0 

M+V.(puu)=Iv.a 


Tt  +  u  .  VP  +  7PV.U  = 


1  v iiiVT)  +  (7 


Re 


where 


is  the  viscous  stress  and 


7-^ooP  =  PT> 

a  =  2p  [i(Vu  +  VuT)  -  i(V  u)l] 


$  =  I(Vu  +  VuT)  :  a 
2 


is  the  viscous  dissipation  function. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


We  assume  a  calorically  perfect  gas  with  the  ratio  of  specific  heats  7(=  Cp/Cv)  =  1.4. 
The  Prandtl  number  Pr  =  0.70,  the  Reynolds  number  Re  is  based  on  8 *,  and  Sutherland’s 
law  is  prescribed  for  the  viscosity  p. 


III.  Primary  Instability  Theory 


In  this  section,  we  consider  the  stability  of  an  infinitesimally  small  primary  disturbance 
superimposed  on  a  steady,  laminar  mean  flow.  The  mean  profiles  are  the  similarity  solution  of 
the  2-d  compressible  boundary  layer  equations  applied  to  a  semi-infinite,  adiabatic  flat  plate. 
They  are  numerically  computed  using  a  spectral  technique  described  in  Streett,  Zang  and 
Hussaini  [15].  The  stability  of  high  Mach  number  flows  is  found  to  be  very  sensitive  to  the 
mean  flow.  This  implies  that  highly  accurate  mean  flow  profiles  are  required  for  meaningful 
stability  calculations  of  high  speed  flows.  Fig.  2a  shows  the  mean  velocity  profiles  versus 
y.  Note  that  in  terms  of  dimensionless  y  the  boundary  layer  thickness  mildly  decreases 
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with  increasing  M^,  although,  of  course,  the  dimensional  boundary  layer  thickness  becomes 
significantly  larger  for  higher  free  stream  Mach  numbers.  Fig.  2b  shows  that  for  MTO  =  4.5, 
there  is  a  substantial  variation  of  the  mean  temperature  profile  across  the  boundary  layer, 
and  that  the  maximum  value  of  the  temperature  profile  is  about  4.5  times  that  of  the  mean 
velocity  profile.  The  marked  change  in  the  temperature  profiles  can  affect  the  stability 
behavior  of  boundary  layers.  The  wall  is  adiabatic,  as  evidenced  by  the  zero  slope  of  the 
temperature  profile. 

Although  the  mean  flow  is  in  general  a  function  of  the  x  and  y  coordinates,  we  invoke, 
in  the  stability  calculations,  the  classical  parallel  flow  approximation  where  the  streamwise 
variation  of  the  mean  flow  is  assumed  negligible  compared  to  its  variation  in  the  normal 
direction.  Consequently  the  primary  disturbance  can  be  written  in  the  modal  form 

q1  =  qi(2/)expW“^-wt)]+c.c.  (7) 

where 

q  =  [u,v,w,T,p],  (8) 

and  the  subscript  1  refers  io  the  primary  wave,  a  and  /3  are  respectively  the  wavenumbers 
in  the  streamwise  and  the  spanwise  directions,  u  is  the  temporal  frequency,  qi(y)  is  the 
complex  amplitude  function  which  determines  the  transverse  structure  of  the  disturbance, 
and  c.c.  denotes  the  complex  conjugate  of  the  preceding  term. 

Two  classess  of  theories  are  prevalent  in  the  study  of  the  stability  characteristics  of  a 
small  disturbance.  In  spatial  theories,  u  is  real  while  a  and  (3  are  complex.  The  real  part 
of  a  is  the  streamwise  wavenumber,  while  the  imaginary  part  of  a  is  the  spatial  growth 
rate.  In  temporal  theories,  a  and  /3  are  real  while  lj  is  complex.  The  real  and  imaginary 
components  of  u  respectively  determine  the  real  frequency  and  the  temporal  growth  rate  of 
the  wave.  Eq.  7  then  defines  a  wave  with  a  real  total  wavenumber  a  =  \Ja2  +  j32  propagating 
with  a  complex  phase  velocity  c  =  w/a  in  the  direction  8  =  tan-1  /3/a.  Temporal  and 
spatial  theories  give  identical  results  on  the  neutral  curve  only  —  where  a,  (3  and  u)  are  all 
purely  real  and  the  growth  rate  is  zero.  When  the  growth  rates  are  sufficiently  small,  an 
approximate  relationship  between  the  temporal  and  the  spa'.ial  grov.  ,h  rate  is  given  by  the 
Gaster  transformation  [5].  In  the  present  ana’ysis,  only  temporally  growing  disturbances  are 
considered. 


Lees  and  Lins  [8]  pioneered  the  theoretical  investigation  of  the  linear  stability  of  in- 
viscid  compressible  flows.  They  found  that  a  sufficient  condition  for  the  existence  of  a 
neutral  subsonic  disturbance  (a  neutral  subsonic  disturbance  decays  monotonically  in  the 
free-stream  with  increasing  y )  is  the  presence  of  an  inflection  point  (which  is  the  point  where 
~  P°  an<^  u°  are  the  mean  density  and  the  mean  velocity).  A  comprehensive 
study  of  the  linear  stability  of  compressible  flows  is  found  in  Mack  [10].  More  recent  re¬ 
sults  on  the  stability  of  compressible  boundary  layers  have  been  obtained  by  Malik  [9].  For 
convenience,  some  of  the  relevant  results  are  restated  here.  By  solving  the  full  disturbance 
equations  (including  the  viscous  and  conductivity  terms)  numerically,  Mack  discovered  that 
several  unstable  modes  may  coexist  in  a  high  Mach  number  compressible  boundary  layer. 
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The  first  mode  is  the  compressible  extension  of  the  TS  mode  found  in  the  incompress¬ 
ible  boundary  layer.  At  low  Mach  numbers,  this  mode  is  primarily  viscous  (the  maximum 
amplification  rate  increases  with  decreasing  Reynolds  number).  It  is  known  that  compress¬ 
ible  boundary  layers  on  insulated  flat  plates  have  a  generalized  inflection  point.  At  higher 
Mach  numbers,  inflectional  instabilities  become  important  because  the  generalized  inflection 
point  is  located  farther  away  from  the  wall.  Therefore,  as  the  Mach  number  increases,  the 
instability  mechanism  progressively  switches  from  the  viscous  to  the  inviscid  type  (where 
amplification  rate  increases  with  increasing  Re),  although  compressibility  has  in  general  a 
stabilizing  influence.  In  fact,  when  M, «,  approaches  about  2.2,  the  viscosity  is  stabilizing  for 
all  Reynolds  numbers,  and  inviscid  instability  modes  are  the  most  unstable.  In  contrast  to 
subsonic  flows,  the  most  unstable  first  mode  is  oblique  for  supersonic  flows. 

Higher  instability  modes  exist  when  the  free  stream  Mach  number  becomes  larger  than 
2.2.  These  modes  do  not  require  a  generalized  inflection  point,  and  they  are  characterized 
by  phase  velocities  close  to  unity  and  high  frequencies  (or  wavenumbers).  The  higher  modes 
depend  on  the  existence  of  a  region  where  the  local  Mach  number  of  the  mean  flow  (relative 
to  the  phase  velocity  c)  is  supersonic.  The  first  of  the  higher  modes  is  called  the  second 
mode.  It  first  appears  when  is  approximately  2.2,  but  is  not  dominant  until  Mm  reaches 
about  3.7.  The  second  mode  reaches  a  peak  amplification  rate  when  M «  4.5.  It  is  the 
most  unstable  of  the  higher  modes,  and  is  most  unstable  when  it  propagates  in  the  mean 
flow  direction. 

The  influence  of  M on  the  eigenfunctions  of  the  first  mode  primary  wave  is  illustrated  in 
Fig.  3  which  shows  that  at  M <*,  =  4.5,  the  maximum  value  of  the  temperature  perturbation 
is  about  15  times  larger  than  the  corresponding  velocity  fluctuation.  This  is  in  contrast  to 
subsonic  flows  where  the  streamwise  velocity  perturbation  is  the  dominant  one.  A  side  by 
side  comparison  between  typical  eigenfunctions  of  the  first  and  second  modes  at  =  4.5 
is  shown  in  Fig.  4.  We  note  that  although  the  vertical  velocity  perturbation  for  the  first 
mode  is  negligible,  it  is  not  insignificant  for  the  second  mode.  Another  feature  of  note  is 
that  the  second  mode  disturbance  decays  much  faster  in  the  freestream  than  does  the  first 
mode  wave. 


IV.  Secondary  Instability  Theory 


The  instability  of  a  base  flow  consisting  of  a  small  amplitude  primary  disturbance  super¬ 
imposed  on  a  parallel  mean  flow  is  called  secondary  instability.  We  begin  with  a  physical 
description  of  how  a  small  monochromatic  primary  disturbance  propagates  downstream  in  a 
boundary  layer.  As  the  disturbance  crosses  the  lower  branch  of  the  neutral  stability  curve, 
its  amplitude  begins  to  grow  exponentially  according  to  linear  (primary)  spatial  instability 
theory.  Assuming  nonlinear  effects  remain  negligible,  the  disturbance  continues  to  amplify 
linearly  until  it  reaches  the  upper  branch  (also  known  as  branch  two)  where  its  local  ampli¬ 
fication  rate  is  again  zero  and  its  amplitude  reaches  a  maximum  value.  Farther  downstream. 
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the  disturbance  begins  to  decay.  When  nonlinear  effects  are  negligible,  it  is  physically  real¬ 
istic  to  use  a  primary  disturbance  which  is  located  near  branch  two  as  input  to  a  secondary 
instability  theory  (since  the  stability  characteristics  of  a  linear  spatial  instability  wave,  which 
achieves  its  maximum  amplitude  at  branch  two,  can  be  accurately  represented  by  a  linear 
temporal  wave  in  the  neighborhood  of  the  neutral  point). 

The  flow  in  the  presence  of  the  primary  disturbance  is  periodic  in  the  streamwise  direction, 
and  hence  its  stability  is  governed  by  Floquet  theory.  For  incompressible  flows,  it  is  known 
that  the  periodic  base  flow  is  very  unstable  to  3-d  secondary  disturbances  [7]. 

We  now  apply  Herbert’s  approach  to  compressible  flows  with  an  oblique  primary  wave. 
When  the  primary  wave  is  oblique,  it  is  appropriate  to  use  a  new-axis,  xx,  which  is  aligned  in 
the  propagation  direction  of  the  primary  wave.  The  Z\  axis  is  perpendicular  to  the  wavevector 
of  the  primary  wave  and  is  in  the  plane  of  the  plate.  Only  when  the  primary  wave  is  two 
dimensional  is  the  fi  axis  oriented  in  the  spanwise  direction 

Specifically,  in  the  (ii,yi  =  y ,  z i)  coordinate  system  (Fig.  1),  the  mean  velocity  profile 
becomes 

u0  =  (u0  cos  0,0,  -uosin0)  (9) 

and  the  primary  velocity  disturbance  takes  the  form 

ui  =  («i>vi,ii>l)e,a(ll_c‘)  +  c.c.  (10) 

where 

iq  =  fix  cos  $  +  uq  sin  0  (11) 

v\  =  hi  (12) 

thi  —  vj\  cos  6  —  hi  sin  9.  (13) 

The  tilde  refers  to  a  quantity  in  the  new  coordinate  system,  the  hat  corresponds  to  a  quantity 
in  the  old  system,  and  c.c.  denotes  the  complex  conjugate  of  the  preceding  term.  Note 
that  there  are  two  effects  associated  with  an  oblique  primary  wave.  First,  there  is  a  mean 
velocity  component  in  the  Z\  direction.  Second,  the  effective  free  stream  Mach  number  in 
the  propagation  direction  of  the  primary  wave  is  reduced  to  M cos  8  from  M^. 

For  a  small  primary  amplitude,  A,  the  base  flow,  say,  for  instance,  the  base  temperature 
Tb,  can  formally  be  expanded  as 

Tb(xuytt)  =  {f0(y)  +  A2T^y)}  +  A{f1(y)  +  A2fl{y)}e'&3*+{A2f}}e2'&i'+0(A4)  (14) 

where  Tq  ,  ,  TJ  are  the  corrections  due  to  the  nonlinear  self-interaction  of  the  primary 

wave.  Similar  expansions  apply  for  other  variables.  When  the  primary  wave  saturates,  its 
amplitude  can  in  principle  be  obtained  from  a  nonlinear  analysis  or  from  direct  numerical 
simulations  [12].  Following  Herbert’s  analysis  of  incompressible  wall  boundary  layers,  we 
assume  that: 
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1.  The  nonlinear  self-interaction  of  the  primary  wave  does  not  significantly  distort  the 
shape  of  the  mean  flow  and  of  the  primary  wave.  For  incompressible  flows,  it  is  found 
that  this  distortion  is  minimal  even  for  strcamwise  velocity  disturbances  of  up  to  10% 
(see  Santos  fl4],  page  11). 

2.  The  amplitude  of  the  primary  wave  remains  constant  during  the  evolution  of  the  sec¬ 
ondary  disturbance.  In  other  words,  the  growth  of  the  primary  wave  is  frozen  by  setting 
the  imaginary  part  of  c,  C{,  to  zero.  The  primary  merely  acts  as  a  catalyst  which  feeds 
energy  from  the  mean  flow  to  the  (secondary)  disturbance.  This  assumption  can  be 
justified  a  posteriori  by  accepting  only  those  secondary  disturbances  which  have  growth 
rates  much  larger  than  that  of  the  primary  wave.  Consequently,  this  assumption  also 
forbids  the  construction  of  neutral  curves  for  the  secondary  disturbance. 

3.  The  amplitude,  £,  of  the  secondary  disturbance  is  assumed  sufficiently  small  to  ensure 
that  linearization  of  the  governing  equations  with  respect  to  z  is  permissible.  A  con¬ 
sequence  of  this  assumption  is  that  the  feedback  from  the  secondary  disturbance  to 
the  primary  wave  is  neglected.  A  complete  analysis  of  this  feedback  mechanism  for 
incompressible  flows  is  detailed  in  Crouch  [2]. 


With  these  assumptions,  the  dependent  variable,  q,  is  decomposed  according  to 

q  =  qb  +  eq2  (15) 

where 

qs,  -  qo(S')  -f  ^{qi(l'K&(il~5rt)  +  c.c}  +  OM2),  (16) 

q0  =  [u0  cos  0,0,  -uosm0,To,  p0]T ,  (17) 

qi  =  [ui,vi,w1,T1,p1]T,  (18) 

Cr  is  the  real  part  of  c,  and  the  subscript  b  indicates  a  base  flow  quantity.  The  viscosity, 
assumed  a  function  of  the  temperature  alone,  is  Taylor  expanded  according  to 

p(T)  =  po  -f  fi'0{ATx  -f  eT2)  +  p.Q{TxT2Ae)  +  0(A2,  e2),  (19) 

where  the  primes  denote  differentiation  with  respect  to  the  mean  temperature. 


Since  qj  is  an  eigenfunction,  the  amplitude  of  the  primary  disturbance  is  not  uniquely 
defined  from  a  linear  stability  analysis.  In  the  (linear)  secondary  instability  analysis  of  El- 
Hady  [6]  and  Nayfeh  [11],  the  eigenfunction  is  normalized  so  that  A  corresponds  directly 
to  the  maximum  r.m.s.  value  of  the  streamwise  velocity  perturbation,  which  is  the  domi¬ 
nant  perturbation  in  subsonic  flows.  For  Mach  numbers  that  are  less  than  2.0,  the  same 
normalization  is  used  here,  i.e., 

max  \iix{y)\2  =  1/2  (20) 

0<y<oo 
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We  caution  the  reader  that  the  above  normalization  may  invalidate  the  use  of  the  shape 
assumption  in  the  supersonic  and  hypersonic  regimes.  For  example,  at  M «,  =  4.5,  the 
maximum  temperature  perturbation  is  an  order  of  magnitude  larger  than  the  corresponding 
maximum  velocity  fluctuation.  This  means  that  if  A  is  defined  in  terms  of  the  velocity 
perturbation,  the  total  magnitude  of  the  temperature  perturbation  AT\  may  no  longer  be 
smaller  than  that  of  the  mean  temperature  To ,  and  consequently  the  neglect  of  terms  of 
0(A2)  and  higher  in  the  construction  of  the  base  flow  (16)  would  not  be  justified.  In  order 
to  avoid  this  possible  inconsistency,  we  define  A,  for  Mach  numbers  >  2,  so  that  it  is  the 
maximum  r.m.s.  value  of  the  largest  perturbation  (which  happens  to  be  the  temperature 
perturbation  in  all  cases  described  here)  i.e., 

max  |qi(t/)|2  =  1/2  (21) 

0<y<oo 

The  above  normalization  ensures  that  when  A  is  small,  the  rms  of  all  primary  variables  are 
also  small.  We  note  in  passing  that  a  normalization  which  is  less  sensitive  to  the  location  of 
the  grid  points  is  by  normalizing  the  primary  eigenfunction  so  that  its  magnitudes  have  a 
specific  area,  for  example, 

r\Uy)\'  =  ll2dy. 

Jo 

We  introduce  a  new  coordinate  system  x  —  x1  —  crt,  y  =  y,z  —  zx  which  moves  with  the 
phase  velocity,  c_,  of  the  primary  wave.  In  this  moving  frame  of  reference,  the  coefficients 
of  the  secondary  disturbance  equations,  which  depend  only  on  the  periodic  base  flow,  are 
independent  of  both  z  and  t.  Thus,  we  can  assume  a  normal  mode  representation  in  z  and 
t  for  the  secondary  disturbance  and  write 


q2  =  q  2(x,y)e"Vfti  (22) 

where  02  is  the  wavenumber  in  the  z  direction.  When  e  is  sufficiently  small,  the  stability 
of  the  secondary  disturbance  is  governed  by  a  system  of  linear  partial  differential  equations. 
These  equations  are  obtained  by  substituting  Eqs.  15-19  into  Eqs.  1-3,  and  keeping  only 
terms  linear  in  e.  Apart  from  a  complicated  y  dependence,  these  equations  are  of  the 
classical  Floquet  type  because  the  coefficients  are  periodic  in  x.  Accordingly,  the  structure 
of  the  secondary  disturbance  in  the  x  direction  has  the  following  general  Floquet  form  (in 
terms  of  a  Fourier  series): 

4  =  ^**  E  (23) 

i=-oo 

where  7  is  the  characteristic  exponent,  e h  is  a  parameter  whose  value  distinguishes  the  various 
types  of  secondary  instabilities,  q2iJ  is  the  shape  function  of  a  mode  whose  wavenumber  in 
the  x  direction  is  given  by  a2  =  e^a/2  +  7  +  ja.  The  fundamental  resonance  is  associated 
with  the  harmonic  series  (e*  =  0)  whose  period  equals  the  primary  wavelength  (i.e.,  27t/^), 
while  the  subharmonic  resonance  is  determined  by  the  subharmonic  series  (e/,  =  1)  whose 
period  is  twice  the  primary  wavelength.  When  0  <  |e^|  <  1,  the  secondary  waves  are  called 
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combination  resonance  modes  (or  detuned  modes).  They  can  have  an  arbitrary  wavenumber 
ip  tne  x  direction,  depending  on  the  value  of  e^. 


We  consider  only  temporal  theory,  which  implies  that  7  —  0  and  that  a  ^  0  is  the  complex 
eigenvalue  to  be  found.  The  real  part  of  o,  aT,  is  the  growth  rate,  and  the  imaginary  part 
of  a,  <7,,  defines  the  frequency  of  the  secondary  disturbance  in  the  moving  coordinate  frame. 
Note  that,  for  a  given  class  of  secondary  waves,  all  the  modes  q2i;  grow  or  decay  at  the  same 
rate  and  travel  with  the  same  frequency  according  to  Eqs.  22-23.  The  growth  rates  are  the 
same  both  in  the  moving  frame  and  in  the  laboratory  frame.  When  a  is  real,  the  secondary 
and  the  primary  disturbances  travel  with  the  same  phase  velocity  in  the  x  direction.  They 
are  phase  locked.  For  complex  cr,  the  phase  velocities  of  the  two  disturbances  are  different. 

Substitution  of  Eqs.  15,  22  and  23  into  the  linearized  disturbance  equations  yields  an 
infinite  set  of  ODE’s  governing  the  stability  ot  the  secondary  disturbance.  Because  the 
governing  equations  for  compressible  flows  have  cubic  convective  nonlinearities  and  because 
the  viscosity  is  a  nonlinear  function  of  the  temperature,  terms  of  0(A2e)  and  higher  are 
present.  These  higher  order  terms  are  neglected  in  the  analysis  because  the  base  flow  is  only 
approximated  to  0(A),  as  are  the  base  flows  of  Herbert  [7],  El  hady  [6],  and  Nayfeh  [11]. 
These  secondary  disturbance  equations  form  a  generalized  complex  eigenvalue  problem 

(axD2  +  A2D1  +  A3j  (f>  —  crA4  4>,  (24) 

where  Dp  indicates  a  p-fold  differentiation  with  respect  to  y.  The  A’s,  which  are  function 
of  the  base  flow  alone,  are  complex,  and  is  the  vector  of  the  unknown  q2  2.  Since  $  is  an 
eigenfunction,  its  amplitude  is  determined  only  up  to  an  arbitrary  complex  constant. 


A  numerical  solution  to  Eq.  24  requires  the  truncation  of  the  Fourier  series  (Eq.  23). 
For  an  N  mode  truncation,  the  size  of  the  A’s  is  5 N  x  5N.  When  N  is  an  odd  integer,  the 
secondary  disturbance  takes  the  form 


q2  =  VW*  £  q2,J(y)e^i. 


N-l 

2 


J=- 


For  an  even  N ,  the  secondary  disturbance  is  represented  by 


(25) 


N-2 

q2  =  ]T  q2j(y)e'J&i  .  (26) 

3=~T 

Accordingly,  to  lowest  order,  the  fundamental  and  the  detuned  disturbances  can  be  repre¬ 
sented  by  3  Fourier  modes  (N  =  3)  according  to 

q2  =  +  q2,0  +  q2Ac'™)  (27) 

and  the  subharmonic  disturbance  by  2  Fourier  modes  (N  =  2)  according  to 

q2  =  e<rte‘^i(q2  _ie-'(«/2)i  +  q20e*(«/2>*)  (28) 
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For  incompressible  flows,  the  aperiodic  q2io  term  in  Eq.  27  is  known  to  give  rise  to  both  a 
mean  flow  distortion  and  a  spanwise  periodic  longitudinal  vortex  structure  associated  with 
the  peak-valley  splitting  in  a  K  type  breakdown  process.  It  turns  out  (see  section  6)  that 
when  the  primary  wave  is  2-d,  the  truncation  levels  defined  by  Eqs.  27-28  are  usually 
sufficient  to  resolve  the  streamwise  structure  of  the  secondary  disturbance,  whether  the  flow 
is  compressible  or  not.  However,  when  the  primary  wavevector  is  inclined  at  a  significant 
angle  to  the  mean  flow  direction,  higher  Fourier  modes  are  required. 

The  stability  calculations  for  compressible  flows  are  substantially  more  computer  inten¬ 
sive  those  of  incompressible  flows  because  all  the  dependent  variables  are  now  coupled  and 
the  number  of  grid  points  required  generally  increases  with  increasing  M A  spectral  repre¬ 
sentation  (see  next  section)  of  Eq.  24  leads  to  a  generalized  matrix  eigenvalue  problem  whose 
eigenvalues  require  0(5  x  N  x  Nc )2  storage  and  0(5  X  IV  x  Nc)2  operations  —  where  Nc  is 
the  number  of  grid  points.  When  the  structure  of  the  compressible  secondary  disturbance  in 
the  x  direction  is  represented  by  more  than  three  Fourier  modes,  it  becomes  impractical  to 
compute  these  eigenvalues  (using  a  global  method)  in  a  parametric  study,  especially  for  the 
higher  Mach  number  flows.  Therefore,  our  computations  are  done  using  Eqs.  27-28,  unless 
otherwise  indicated.  The  effect  of  the  neglected  higher  Fourier  modes  on  a  is  only  examined 
in  some  selected  cases,  although  some  insight  into  the  convergence  of  the  Fourier  series  can 
also  be  obtained  from  DNS. 

The  appropriate  boundary  conditions  for  the  secondary  disturbances  are 

(u2,  v2,w2,T2)  =  0  at  y  =  0  (29) 

(u2,v2,w2,T2)  -»  0  as  y  ->  oo  (30) 

The  boundary  conditions  on  the  density  disturbances  are  not  explicitly  imposed  since  they 
satisfy  the  continuity  equation.  We  next  outline  the  numerical  algorithm. 


V.  Numerical  Method 


The  system  of  homogeneous  equations  Eq.  24  plus  the  homogeneous  boundary  conditions, 
Eqs.  29-30,  lead  to  a  dispersion  relation  of  the  form 

a  =  a(a2, /?2)-  (31) 

The  eigenvalue,  a ,  is  obtained  numerically  using  a  spectral  collocation  technique  based  on 
Chebyshev  polynomials.  The  collocation  points  for  all  the  dependent  variables  are  the  Gauss- 
Lobatto  points.  These  points,  which  are  the  extrema  of  the  highest  Chebyshev  polynomial 
retained,  are  defined  by 

Vj  -  cos  ^  j  =  0,1,...,  W  (32) 

where  Af  is  highest  order  Chebyshev  polynomial.  In  order  to  provide  an  adequate  resolution 
near  the  wall  and  the  critical  layer,  the  physical  domain  y  is  mapped  onto  the  Chebyshev 
domain  tj  through  two  different  transformations. 
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The  first  transformation,  upon  which  most  of  the  results  reported  in  this  paper  are 
based,  maps  the  finite  physical  domain  y  £  [0,  ymax]  onto  rj  £  [—1, 1]  through  a  combination 
of  hyperbolic  tangent  and  algebraic  stretchings  according  to 


v- Vo 
At] 


(33) 


where 


j/l/23/max(l  +  VO 
1/max  ^(l/max  2  J/i  y2  ) 

dip 

=  dy 'u=VoAy°' 


(34) 


(35) 


Here  te  is  a  concentration  parameter  which  clusters  nodes  about  ipQ.  If  tt  =  0,  there  is 
no  stretching  due  to  the  influence  of  the  hyperbolic  tangent  term.  Ayo  is  the  width  of 
the  concentration  region  in  the  physical  domain.  The  far  field  boundary  of  the  discretized 
physical  domain  is  located  at  ymax,  and  half  of  the  nodes  are  distributed  between  y  =  0 
and  y  =  yi/2-  Dirichlet  conditions  are  imposed  on  both  the  primary  and  the  secondary 
disturbances  at  y  =  ymax  where  all  variables,  except  the  density,  are  set  to  zero.  The 
determination  of  the  adjustable  constants  ymax,  yi/2,  tt  and  A y0  is  empirical;  more  details 
are  given  in  [4].  For  <  1.6,  ymax  =  100,  y1/2  =  2,  and  te  =  0.  When  Mm  =  4.5,  two 
sets  of  parameters  are  used;  for  first  mode  calculations,  runs  are  performed  with  ymax  — 
30,  yx/2  =  1,  t€  =  0.8,  y0  =  1.2  and  Ay0  =  0.4,  while  ymax  =  15,  y1/2  =  1,  tt  =  0.8,  y0  =  1.2 
and  Ayo  —  0.4  are  for  second  mode  calculations.  A  lower  value  of  ymax  for  second  mode 
calculations  reflects  the  faster  decay  of  the  second  mode  eigenfunctions  in  the  free  stream 
compared  to  that  of  the  first  mode  waves. 


The  second  mapping  transforms  y  £  [0,oo]  onto  r)  £  [1,0]  via 


V  = 


h  .  1  y\ 
2bTT)  +  2exp[-kl 


(36) 


When  n  =  1  and  ij  =  0,  the  mapping  reduces  to  the  classical  exponential  mapping.  For 
given  lx  and  l2,  n  controls  the  number  of  grid  points  in  the  neighborhood  of  the  critical  point 
(near  y  =  1).  Increasing  n  adds  more  nodal  points  around  the  critical  layer  region.  The 
parameters  l2  and  n  redistribute  some  of  the  nodal  points  near  y  =  0  to  the  critical  layer. 
The  point  y  =  0  is  mapped  onto  77  =  1,  while  the  point  y  =  00  is  mapped  onto  77  =  0  (not 
a  collocation  point).  All  variables  are  expanded  in  odd  polynomials  to  implicitly  satisfy  the 
boundary  conditions  at  y  =  00.  Calculations  are  performed  with  l\  =  7,  l2  =  7,  and  n  =  3. 

A  global  method,  based  on  a  modified  QZ  algorithm  [16],  is  used  to  obtain  all  the  eigen¬ 
values  of  the  discrete  system  of  linear  equations.  We  then  select  the  most  unstable  eigenvalue 
and  use  it  as  input  to  a  local  procedure,  which  is  based  on  an  inverse  Rayleigh  solver.  The 
local  method  generates  a  more  accurate  eigenvalue  and  its  corresponding  eigenfunction.  The 
eigenvalues  are  at  least  accurate  to  5  decimal  places  for  primary  eigenfunctions  and  to  4  dec¬ 
imal  places  for  a  secondary  eigenfunctions.  Since  the  accuracy  of  the  secondary  disturbance 


11 


depends  on  the  accuracy  of  the  primary  wave  but  not  conversely,  a  more  stringent  conver¬ 
gence  criterion  is  imposed  on  the  primary  wave  to  ensure  that  its  derivatives  are  determined 
with  sufficient  accuracy.  In  all  cases,  the  ratio  of  the  maximum  of  the  absolute  value  of  the 
last  three  Chebyshev  coefficients  over  the  maximum  of  the  absolute  value  of  the  first  half 
Chebyshev  coefficients  is  always  less  than  10'4;  this  procedure  should  exclude  any  spurious 
modes,  unless  they  happen  (very  unlikely)  to  be  well  resolved. 
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VI.  Results  and  Comparison  with  Direct  Numerical 

Simulations 


A.  Code  validation 

In  the  absence  of  experimental  data  for  compressible  flows,  the  secondary  instability  code  is 
validated  against  the  established  incompressible  results  of  Herbert  et  al  [7]  and  against  re¬ 
sults  from  direct  numerical  simulations  of  the  full  compressible  Navier-Stokes  equations. 
For  convenience,  the  verification  for  incompressible  flow  is  reported  using  the  reference 
length  L  =  (1/^,1* /u^,)1/2  instead  of  8*.  The  parameters  used  (in  Herbert’s  units)  are 
a  =  0.20335, =  0,  A  =  1  %,A>  =  0.2,  =  10"6,  Re  =  606  and  Too  =  520 °R.  The  results 

are  tabulated  below: 


Detuning 

Number  of 
Fourier  Modes 

45  points 

55  points 

KEQ 

123 

mm 

■m 

eh  =  1 

2 

8.1177 

0 

8.1177 

0 

subharmonic 

3 

8.1568 

0.0145 

8.1568 

0.0144 

mode 

4 

8.1958 

0 

8.1958 

0 

eh  =  0.5 

2 

6.1811 

1.6457 

6.1811 

1.6456 

detuned 

3 

6.4067 

1.4818 

6.4067 

1.4816 

mode 

4 

6.4185 

1.4711 

6.4183 

1.4710 

c*  =  0 

2 

-1.7356 

2.0434 

fundamental 

3 

3.6388 

0 

3.6387 

0 

mode 

4 

3.6498 

-0.0102 

3.6495 

-0.0103 

Table  1 

A  side  by  side  comparison  between  the  above  table  and  Table  1  of  Herbert  et.  al.  [7]  shows 
that  the  eigenvalues  are  in  agreement  to  at  least  5  decimal  places.  Except  for  Table  1,  all 
computations  are  based  on  6*. 

The  verification  for  compressible  flows  is  performed  using  a  code  described  in  detail  in 
Erlebacher  and  Hussaini  [4].  For  completeness,  a  brief  description  of  the  method  is  outlined 
here.  Periodic  boundary  conditions  in  x  and  z  are  imposed  where  x  is  in  the  propagation  di¬ 
rection  of  the  2-d  or  3-d  primary  wave.  The  computation  domain  has  dimensions  L*  =  27 t/q* 
and  L-z  =  2ir/(3  in  the  x  and  z  directions  respectively,  where  a*  =  a  for  the  fundamental 
instability  and  a*  =  a/ 2  for  the  subharmonic  instability.  In  the  normal  direction,  the  finite 
physical  domain  y  6  [0,ym<Jx]  is  mapped  onto  rj  (E  [—1,1]  according  to  Eqs.  33-35.  The 
spatial  discretization  is  based  on  a  spectral  collocation  method  with  Fourier  expansions  in 
the  x  and  z  directions,  and  Chebyshev  expansions  in  the  normal  direction,  y  =  y.  Velocities 
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satisfy  no-slip  boundary  conditions  and  the  wall  is  adiabatic.  Dirichlet  boundary  conditions 
are  imposed  at  y  —  ymax .  Initially,  the  flow  consists  of  a  parallel  mean  flow  plus  the  primary 
and  the  secondary  disturbances.  Time  marching  is  achieved  with  a  third  order  Runge  Kutta 
method.  The  growth  rate  gji  of  the  disturbances  is  monitored  through 

1  dEji 
9,1  =  2  Eji  dt 


where 


Eji{t)  =  /  (M2  +  bj.il2  +  \wjA2)dy, 

Jo 


(38) 


Uji,  Vji  and  Wji  are  the  Fourier  coefficients  of  the  u,  v  and  w  velocity  components,  respectively. 
The  subscripts  j  and  l  denote  the  Fourier  components  in  the  x  and  z  direction,  respectively. 


Comparisons  between  theory  and  DNS  are  performed  at  Mach  1.6  (Fig.  5a)  and  Mach  4.5 
(Fig.  5b).  At  Moo  =  1.6,  the  parameters  are  Re  =  1675,  Too  =  520°72,  (a,  8)  =  (0.24,10°), 
f 32  =  0.4,  e/,  =  0,  N  =  3,  A  =  0.015,  and  e  =  0.0015.  The  calculations  at  JlfM  =  4.5  are  for 
Re  =  10000,  and  T^  =  110°R,  (a,  8)  =  (2.52,0°),  (S2  =  2.1,  eh  =  1,  N  =  2,  A  =  0.06,  and 
e  =  0.006.  Here  t  is  the  initial  amplitude  of  the  secondary  disturbance  used  in  the  DNS,  and 
time  is  normalized  with  respect  to  the  period  of  the  primary  wave.  The  following  two  tables 
represent  respectively  samples  of  each  comparison  at  Mw  =  1.6  (Table  2)  and  =  4.5 
(Table  3): 


time 

primary  growth  rate  xlO4 

secondary  growth  rate  x  102 

theory 

DNS 

theory 

DNS 

0.000 

1.919 

1.911 

1.127 

1.112 

0.010 

1.919 

1.799 

1.127 

1.111 

0.023 

1.919 

1.947 

1.127 

1.112 

Table  2 


time 

primary  growth  rate  xlO3 

secondary  growth  rate  xlO2 

theory 

DNS 

theory 

DNS 

0.000 

2.867 

3.023 

2.506 

2.500 
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2.506 
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Fig.  5a  represents  the  dominant  fundamental  mode  at  =  1.6  and  A  =  1.5%,  which 
arises  from  a  primary  disturbance  propagates  at  10°  to  the  mean  flow  direction,  while  Fig.  5b 
represents  the  most  amplified  subharmonic  based  on  a  2-d,  second  mode  primary  disturbance 
with  A  —  6%.  The  figures  show  that  the  theoretical  growth  rates  of  primary  and  secondary 
disturbances  compare  well  with  the  growth  rates  obtained  by  DNS.  There  is  a  small  difference 
between  the  Prandtl  number  used  in  theory  (0.70)  and  in  DNS  (0.72).  This  difference, 
although  small  and  unintentional,  may  partly  account  for  the  small  discrepancies  between 
theoretical  and  numerical  results.  For  both  the  Mach  1.6  and  Mach  4.5  cases,  the  secondary 
disturbances  are  stable  in  the  absence  of  the  primary  disturbance.  However,  the  presence 
of  the  finite  amplitude  primary  disturbance  triggers  a  strong  secondary  instability  with  the 
growth  rate  of  the  secondary  disturbance  far  exceeding  that  of  the  primary. 

B.  Parametric  studies 

The  free  stream  Mach  number  M ^  affects  secondary  instability  directly  through  the  (sec¬ 
ondary)  disturbance  equations  and  indirectly  through  modification  of  the  mean  flow  and 
primary  wave.  The  influence  of  Af^  on  secondary  instability  is  studied  in  this  section. 


VI.2.1.  Subsonic  flows 

We  begin  with  a  study  of  the  secondary  instability  at  M ^  =  0.8.  Fig.  6  depicts  the  growth 
rate  versus  the  spanwise  wavenumber  for  the  three  types  of  secondary  instability  waves  at 
Re  =  1675,  Too  =  520°/2,  a  =  0.29,  (3  =  0,  and  two  primary  wave  amplitudes:  A  -  1.5% 
and  A  =  0.5%.  The  parameters,  Re  and  Too,  are  approximately  the  same  as  those  given  in 
[6].  The  primary  wave  is  2-d,  slightly  unstable,  and  is  located  near  branch  two  of  a  neutral 
curve  with  a  growth  rate  of  0.0006.  Although  the  subharmonic  and  the  fundamental  dis¬ 
turbances  travel  with  the  phase  velocity  of  the  primary,  the  detuned  disturbance  does  not. 
Additionally,  the  following  salient  features  may  be  extracted  from  this  figure:  (1)  There  is  a 
wide  band  of  highly  unstable  3-d  secondary  disturbances  whose  growth  rates  increase  with 
increasing  primary  wave  amplitude.  As  /32  — >  0,  the  secondary  disturbances  rapidly  become 
stable.  Since  2-d  (i.e.,  =  0)  secondary  disturbances  are  unstable  in  free  shear  layers  [13], 

we  conclude  that  the  presence  of  the  wall  at  y  =  0  can  damp  the  secondary  disturbances  with 
long  spanwise  wavelengths.  (2)  There  is  a  preferred  band  of  spanwise  wavenumbers  within 
which  the  secondary  growth  rate  has  a  local  maximum.  Increasing  the  primary  wave  ampli¬ 
tude  mildly  increases  the  preferred  spanwise  wavenumber  and  significantly  widens  the  range 
of  /?2  over  which  a  secondary  disturbance  is  unstable.  At  A  =  1.5%,  the  spanwise  wavenum¬ 
ber  of  the  most  amplified  subharmonic  approximately  equals  the  streamwise  wavelength  of 
the  primary  wave.  As  the  type  of  instability  changes  from  subharmonic  to  fundamental, 
(i.e.,  as  decreases  from  1  to  0),  the  preferred  spanwise  wavenumber  increases  slightly.  (3) 
The  subharmonic  modes  are  the  most  unstable,  followed  by  the  detuned  modes  and  then 
the  fundamental  modes.  The  preceding  result  is  consistent  with  the  classical  scenario  of  th< 
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H-type  breakdown  in  which  the  subharmonic  modes  dominate  when  the  primary  amplitude 
is  low. 

The  subharmonic  growth  rates  from  Fig.  6  are  then  compared  to  two  sets  of  incompress¬ 
ible  results.  Fig.  7a  illustrates  the  first  comparison  using  Re  =  1675,  T =  520°R,  a  =  0.29, 
and  0  =  0.  The  primary  growth  rates  at  Moo  =  0.0  and  M „  =  0.8  are  0.0028  and  0.0006 
respectively.  For  clarity,  let  0max  be  the  spanwise  wavenumber  of  the  most  amplified  subhar¬ 
monic.  An  interesting  feature  is  that  subsonic  compressibility  destabilizes  the  subharmonic 
waves  with  02  >  0max.  At  A  =  0.5%,  compressibility  appears  to  have  a  stabilizing  influence 
for  the  disturbances  with  02  <  fimax- 

The  second  set  of  calculations  is  performed  using  a  =  0.29  and  two  different  Reynolds 
numbers:  Re  =  2700  for  =  0.0  and  Re  =  1675  for  M^  =  0.8.  These  Reynolds  numbers 
are  chosen  to  ensure  that  the  primary  growth  rates  at  both  Mach  numbers  are  approximately 
the  same  (cjj  =  0.0006).  The  results  for  the  subharmonic  growth  rate  versus  the  spanwise 
wavenumber  are  plotted  in  Fig.  7b.  In  contrast  to  the  results  in  Fig.  7a,  compressibility  now 
has  a  stabilizing  influence.  Therefore,  although  there  is  a  strong  growth  of  3-d  secondary 
disturbances  at  subsonic  Mach  numbers,  the  particular  effect  of  compressibility  on  local 
subharmonic  growth  rates  is  unclear  (we  remark  that  the  growth  rates  have  been  scaled  with 
ulo/S*  which  changes  with  Ma Q).  It  can  either  have  a  stabilizing  or  a  destabilizing  influence, 
depending  on  how  the  comparison  with  incompressible  flow  is  made.  Such  a  comparison  is 
only  meaningful  when  the  following  conditions  are  met  (i)  when  all  quantities  are  expressed 
in  dimensional  form,  (ii)  when  the  free  stream  disturbances  are  coupled  with  the  primary 
disturbances  (i.e  a  receptivity  problem),  and  (iii)  when  the  total  growth  of  the  mean  flow 
and  the  primary  disturbance  is  incorporated  into  the  calculations.  The  point  is  that  since 
the  primary  wave  amplitude  (which  is  responsible  for  the  onset  and  growth  of  the  secondary 
disturbance  and  which  is  in  itself  sensitive  to  the  free-stream  enviroment)  cannot  be  uniquely 
determined  from  stability  analysis,  it  is  meaningless  to  make  any  general  conclusions  on  the 
overall  influence  of  compressibility  on  secondary  instabilities.  Henceforth  conclusions  about 
the  secondary  instabilities  in  this  paper  are  local  in  nature  and  are  based  on  the  primary 
wave  amplitudes  (and  other  flow  parameters)  which  are  assumed  known. 

VI.2.2.  Supersonic  flows 

The  subharmonic  and  the  fundamental  growth  rates  at  Ma 0  =  1.6,  Re  =  1675,  T «,  =  520°  R, 
a  =  0.24,  and  0  =  0  are  plotted  against  02  in  Fig.  8  and  Fig.  9  respectively.  The  primary 
mode  is  slightly  unstable,  located  near  branch  two  of  the  neutral  curve.  There  is  a  broad¬ 
band  of  3-d  secondary  disturbances  whose  growth  rates  are  large  (compared  to  the  maximum 
primary  growth  rate  which  approximately  equals  0.003)  and  increase  with  increasing  A. 
The  results  lead  us  to  believe  that  secondary  instability  mechanisms  at  Mw  =  1.6  play  a 
significant  role  in  the  boundary-layer  transition  and  that  the  process  of  laminar  breakdown 
is  similar  to  the  mechanisms  found  in  incompressible  flows.  The  spanwise  wavelength  of 
the  most  unstable  secondary  disturbance  is  about  60%  of  the  streamwise  wavelength  of  the 


primary  disturbance,  and  this  wavelength  is  almost  insensitive  to  the  primary  amplitude. 

Fig.  10  shows  the  influence  of  the  primary  amplitude  A  on  aT.  The  parameters  used  are 
identical  to  those  of  Figs.  8  and  9  when  /?2=  0.4  (i.e.,  the  spanwise  wavenumber  of  the  most 
dangerous  secondary  disturbance).  Here  the  fundamental  instability  is  stronger  than  the 
subharmonic  instability  when  A  is  larger  than  2.5%,  while  the  converse  is  true  for  smaller 
values  of  A. 

Next  the  secondary  instability  is  studied  M =  4.5 ,  T,*,  =  110°R,Re  =  10000,  and  two 
values  of  a  :  0.6  for  a  first  mode  calculation  and  2.52  for  a  second  mode  calculation  for  four 
primary  amplitudes  (A  =  0%,  2%,  3%,  6%).  The  growth  rates  of  the  2-d  primary  waves  for 
the  first  mode  and  the  second  mode  calculations  are  0.002087  and  0.003185,  respectively. 
The  first  mode  primary  wave,  which  is  the  same  as  that  reported  in  [3],  is  the  most  amplified 
2-d  mode.  The  second  mode  primary  wave  is  chosen  so  that  it  is  located  near  branch  two. 
For  these  parameters,  no  evidence  of  the  fundamental  secondary  instability  can  be  found. 
The  subharmonic  instability  results  are  shown  in  fig.  11a  (first  mode)  and  fig.  lib  (second 
mode). 

In  Fig.  11a,  it  is  interesting  to  note  that  when  the  primary  amplitude  is  non-zero,  there 
are  two  humps  in  the  curve  of  ar  versus  /32-  Hump  one  has  a  complex  a  whereas  hump 
two  has  a  real  a.  The  growth  rates  of  the  unstable  modes  whose  spanwise  wavelengths  are 
longer  than  the  streamwise  wavelength  of  the  primary  mode  are  marginally  affected  by  the 
primary  amplitude;  in  other  words,  these  waves  are  essentially  the  3-d  primary  (subharmonic) 
instability  modes.  The  unstable  waves  on  hump  two  are  secondary  instability  waves.  While 
increasing  the  primary  amplitude  marginally  stabilizes  hump  one,  it  significantly  destabilizes 
hump  two.  At  A  =  6%,  the  growth  rates  of  the  most  amplified  primary  (subharmonic)  and 
secondary  (subharmonic)  modes  are  comparable. 

Fig.  lib  demonstrates  that  the  subharmonic  instabilities  which  originate  from  the  sec¬ 
ond  mode  primary  wave  differ  substantially  from  the  corresponding  first  mode  case.  Here 
<jT  versus  f32  curve  has  only  one  maximum.  The  eigenvalue  a  is  real  for  all  values  of  /32, 
and  the  subharmonic  in  the  absence  of  the  primary  disturbance  is  stable.  Increasing  the 
primary  wave  amplitude  destabilizes  the  subharmonic  and  dramatically  widens  the  range 
of  spanwise  wavenumber  over  which  a  subharmonic  is  unstable.  The  most  amplified  sub¬ 
harmonic,  whose  spanwise  wavenumber  @2  ~  2.1  is  rather  insensitive  to  A,  has  a  spanwise 
wavelength  approximately  equal  to  80%  of  the  streamwise  wavelength  of  the  primary.  As  A 
changes  from  2%  to  6%,  the  growth  rate  of  the  most  amplified  subharmonic  is  increased  by 
135%  (from  0.0075  to  0.0251).  Since  the  maximum  growth  rates  of  both  the  primary  and 
the  secondary  disturbances  for  the  second  mode  calculations  are  considerably  higher  than 
the  corresponding  first  mode  calculations,  we  conjecture  that  boundary-layer  transition  at 
Moo  =  4.5  will  occur  (in  all  likelihood)  via  second  mode  disturbances. 

The  convergence  of  the  Fourier  modes  for  the  most  unstable  subharmonic  mode  from 
Fig.  lib,  together  with  the  fundamental  and  the  detuned  modes,  is  performed.  (To  the 
authors  best  knowledge,  such  a  convergence  study  has  never  been  done  for  compressible 
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secondary  instability  theory).  As  evidenced  from  Table  4,  the  stream  wise  structure  of  the 
subharmonic  disturbance  can  be  accurately  captured  by  using  only  two  Fourier  modes,  just 
as  in  incompressible  flows. 
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Fig.  12  shows  <7r  as  a  function  of  the  detuning  parameter  at  =  4.5  using  3  Fourier 
modes  (i.e.,  N  =  3),  a  =  2.52,  0  =  0,  A  =  0.06,  =  110°fl,  Re  =  10000,  and  /32  =  2.1. 

It  reveals  that  aT  is  a  continuous  function  of  e*  and  that  there  is  a  broad  band  of  highly 
unstable  modes.  The  subharmonic  mode  has  the  largest  growth  rate,  while  the  fundamental 
mode  is  stable. 

The  corresponding  influence  of  the  detuning  parameter  on  the  frequency  shift  <T;  is  illus¬ 
trated  in  Fig.  13.  While  the  subharmonic  mode  has  real  a  (which  implies  that  it  travels 
synchronously  with  the  primary  mode),  the  detuned  and  fundamental  modes  have  complex 
a. 


Fig.  14  shows  the  influence  of  the  primary  amplitude  A  on  the  subharmonic  growth  rate. 
The  parameters  involved  are  identical  to  those  of  Fig.  lib  with  /?2  =  2.1.  We  note  that  <jt 
increases  with  A,  and  that  the  subharmonic  modes  are  unstable  only  when  A  is  larger  than 
1.2% 


VI.2.3.  Secondary  Eigenfunctions 

Our  previous  results  indicate  that  the  primary  wave  can  trigger  a  strong  growth  of  secondary 
disturbances.  We  proceed  next  to  examine  the  influence  of  compressibility  on  the  secondary 
eigenfunctions,  which  determine  the  structure  of  the  secondary  disturbance  normal  to  the 
plate. 
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The  fundamental  eigenfunctions  at  M =  1.6  associated  with  wavenumbers  a2  =  0  and 
a2  —  0-24  are  displayed  in  Figs.  15a  and  15b  respectively.  These  eigenfunctions  correspond  to 
the  dominant  fundamental  instability  with  a  2-d  primary  mode  at  =  1.6  and  A  =  1.5%. 
The  results  are  obtained  using  Re  =  1675,  Too  =  520°  R,  a  =  0.24,  (3  =  0,  /32  =  0.4,  45 
grid  points  and  3  Fourier  modes.  The  two  most  dominant  components  of  the  secondary 
eigenfunctions,  in  order  of  importance,  are  the  streamwise  velocity  components  u2  0  and  u2  j. 
Both  components  peak  near  the  critical  point  (i.e.,  y  ~  1.0).  In  contrast,  at  this  Mach 
number,  the  strongest  primary  eigenfunction  is  the  temperature  component. 

Fig.  15a  shows  that  the  ^-distributions  of  the  temperature,  density  and  spanwise  velocity 
eigenfunctions  all  have  two  local  maxima  across  the  boundary  layer.  4  comparison  with  the 
corresponding  incompressible  case  in  Fig.  16a  reveals  two  compressibility  effects.  First,  at 
higher  Moo  the  peaks  of  the  eigenfunctions  are  pushed  away  (in  dimensional  units)  from 
the  wall.  For  example,  while  u2)0  peaks  near  y *  =  0.76*  (=  t/q)  at  —  0,  it  peaks 
near  y*  —  0.96*  (=  y{)  at  M*  =  1.6.  Since  6*  =  2.72^/i/^i* /u*^  at  Mm  =  1.6  and 

6*  =  1.72^  x*/ula  at  Moo  =  0,  it  is  clear  that  yf  is  about  twice  as  large  as  y£  —  recall  that 
the.  superscript  *  indicates  a  dimensional  quantity.  The  second  effect  of  compressibility  is  to 
make  the  spanwise  velocity  eigenfunction  w2<0  become  insignificant  relative  to  the  streamwise 
velocity  eigenfunction  u2io- 

Fig.  15b  shows  that  the  eigenfunctions  (except  for  u2>i)  associated  with  a2  =  0.24  at 
Moo  =  1.6  are  characterized  by  a  major  peak  across  the  boundary  layer.  The  |u2il  |  component 
has  two  peaks:  one  very  small  one  at  y  «  0.01  (not  clearly  shown  in  the  figure)  and  the 
other  at  y  «  0.9. 

Fig.  16b  shows  that  although  the  velocity  eigenfunctions  at  Af*  =  0  are  qualitatively 
similar  to  Fig.  15b,  there  are  some  noticeable  differences.  For  example,  the  ratio  of  the 
maximum  value  of  ]iy2il|  to  the  maximum  value  of  |u2,i|  is  about  0.8  at  =  0,  whereas 
this  ratio  becomes  0.5  at  =  1.6. 

The  influence  of  compressibility  on  the  subharmonic  eigenfunctions  is  examined  at  Mach 
numbers  0,  1.6  and  4.5,  and  is  illustrated  in  Figs.  17-19.  These  figures  show  that  compress¬ 
ibility  introduces  temperature  and  density  fluctuations  which  become  increasingly  significant 
at  higher  Mach  numbers.  Since  6*(=  10.5yt'*x*/u£o  at  =  4.5)  increases  rapidly  with 
increasing  M0 <,,  the  figures  confirm  that  the  vertical  location  (in  dimensional  unit)  of  the 
subharmonic  perturbations  peaks  are  pushed  towards  the  free  stream  by  increasing  Moo- 
Moreover,  the  eigenfunctions  at  =  4.5  (centered  near  the  critical  point)  are  confined 
to  a  much  narrower  region  than  those  at  lower  Mach  numbers.  Furthermore,  although  the 
streamwise  velocity  perturbation  is  the  dominant  perturbation  at  =  0  and  M^  =  1.6, 
the  temperature  perturbation  becomes  the  largest  one  at  M0 0  =  4.5. 
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VI.2.4.  Reynolds  number  effect 


We  determine  here  the  Reynolds  number  beyond  which  the  inviscid  secondary  disturbance 
equations  can  capture  the  correct  instability  characteristics  of  secondary  disturbances.  To 
this  effect,  we  vary  the  Reynolds  number  (only)  in  the  secondary  disturbance  equations.  A 
calculation  is  performed  for  the  most  amplified  secondary  disturbance  at  M «,  =  4.5  and 
A  =  6%.  The  parameters  are  T «,  =  110°R,  a  =  2.52 ,/3  =  0,  e/,  =  1,  N  =  2,  and  02  =  2.1. 
The  Reynolds  number  used  for  the  primary  mode  is  1  x  104,  while  it  ranges  from  1  x  104  to 
1.49  xlO7  for  the  secondary  mode.  The  results  are  shown  in  Fig.  20.  We  note  that  the  second 
mode  subharmonic  is  inviscid  in  character  in  the  sense  that  aT  increases  (monotonically)  with 
Re  and  asymptotes  to  an  inviscid  limit.  The  figure  demonstrates  that  inviscid  results  can  be 
obtained  from  the  viscous  problem  in  the  limit  of  vanishing  viscosity.  However,  the  growth 
rate  of  the  subharmonic  can  be  obtained  within  1%  accuracy  from  inviscid  calculations  only 
for  Reynolds  numbers  higher  than  a  million.  For  example,  ar  at  Re  =  1  x  104  is  only  about 
45%  of  the  inviscid  growth  rate. 

The  corresponding  effects  of  Reynolds  number  on  the  subharmonic  eigenfunctions  are 
depicted  in  Figs.  19,  21  and  22.  While  the  viscosity  minimally  affects  the  thermodynamic 
variables  t2<\  and  p2,i  (which  are  the  dominant  perturbations),  it  plays  an  appreciable  dif¬ 
fusive  role  on  the  streamwise  velocity  eigenfunction  (u2il)  and  the  the  spanwise  velocity 
eigenfunction  (u>2il). 


VI.2.5.  Primary  wave  angle  effect 

Recall  that  the  most  unstable  primary  disturbance  in  supersonic  and  hypersonic  flows  can  be 
oblique  (with  respect  to  the  mean  flow  direction).  The  influence  of  the  primary  wave  angle 
9  on  <7r  is  examined  at  =  1.6,  Re  =  1675,  T ^  =  520°R,  a  =  0.24,  and  A  =  1.5%.  First 
small  to  moderate  values  of  9  are  examined  for  9  =  0°,  5°,  10°,  15u.  The  primary  modes  are 
slightly  unstable  with  growth  rates  varying  from  3.7  x  10-5  at  8  =  0  to  1.42  x  10~3  at  9  —  45°. 
Since  the  secondary  instabilities  of  a  2-d  primary  wave  can  be  accurately  determined  using 
the  lowest  order  of  approximation  defined  by  Eqs.  27-28,  we  assume  that  this  approximation 
will  not  incur  significant  errors  in  the  determination  of  cr  for  the  aforementioned  values  of  9 ; 
consequently,  the  convergence  of  a  in  terms  of  the  Fourier  modes  is  not  verified.  In  Figs.  23- 
24,  aT  is  plotted  against  /32.  It  is  found  that  when  9  ^  0°,  the  secondary  disturbances  do  not 
travel  with  the  phase  velocity  of  the  primary  disturbance.  The  figures  show  that  while  the 
strongest  subharmonic  mode  occurs  when  the  primary  disturbance  propagates  in  the  mean 
flow  direction,  the  most  unstable  fundamental  mode  occurs  when  9  m  10°.  When  9  >  10°, 
there  are  two  local  maxima  in  the  curve  of  aT  versus  /?2.  Because  the  second  maximum  is 
relatively  smaller  than  the  first  maximum,  it  will  probably  have  no  significant  impact  on 
boundary-layer  transition.  It  is  also  worth  noting  that  the  growth  rates  of  the  unstable 
waves  with  small  @2  are  only  very  weakly  dependent  on  the  orientation  of  the  primary  wave, 
especially  for  small  values  of  9.  Since  neither  the  largest  primary  growth  rate  nor  the  largest 
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fundamental  growth  rate  corresponds  to  a  2-d  primary,  it  is  inappropriate  to  examine  the 
fundamental  instability  based  on  a  2-d  primary  wave,  at  least  at  Moo  —  1-6. 


It  is  known  that  at  M <*,  =  1.6,  the  most  unstable  primary  wave  occurs  at  6  =  45°.  Since 
this  value  of  6  is  appreciably  high,  it  is  appropriate  to  first  perform  convergence  studies  of 
the  effect  of  the  higher  fourier  modes  on  secondary  instabilities.  The  results,  together  with 
the  corresponding  case  for  9  =  0°,  are  illustrated  in  Table  5. 
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We  see  that  when  the  primary  is  2-d,  the  first  two  Fourier  modes  can  capture  a  correctly  to 
3  decimal  places.  However,  when  6  =  45°,  the  first  two  Fourier  modes  are  not  sufficient  to 
accurately  determine  a.  The  first  3  Fourier  modes  determine  a  correctly  to  2  decimal  places. 
The  results  from  Table  5  can  easily  be  explained  by  the  corresponding  eigenfunctions.  The 
normalized  streamwise  velocity  eigenfunctions,  which  are  the  most  dominant  components, 
are  shown  in  Fig.  25  for  0  =  0°  and  Fig.  26  for  9  —  45°.  When  9  =  0°,  the  eigenfunctions 
associated  with  a 2  =  ±a/2  are  the  most  important  components,  while  those  associated  with 
higher  values  of  a2  are  increasingly  negligible.  Consequently,  a  can  be  accurately  determined 
by  using  only  the  first  two  Fourier  modes.  When  9  =  45°,  the  eigenfunctions  associated  with 
a2  =  a/2  and  a2  =  3a/2  are  respectively  the  two  most  important  components.  Other 
eigenfunctions  are  relatively  insignificant.  Therefore,  a  can  be  determined  with  sufficient 
accuracy  using  the  modes  a2  =  a/ 2  and  a2  —  3a/2,  but  not  using  the  components  with 
a2  =  ±a/2.  Finally,  it  is  interesting  to  note  that  the  eigenfunctions  with  a2  =  — 3d/2  is 
practically  negligible.  This  explains  why  a  does  not  change  when  N  is  increased  from  N  =  3 
to  IV  =  4  in  Table  5.  Hence,  caution  must  be  taken  when  performing  the  convergence  studies 
on  the  eigenvalue  a. 
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VII.  Conclusions 


A  fully  spectral  code  has  been  developed  to  study  the  linear  secondary  instabilities  of 
compressible  boundary  layers  on  insulated  flat  plates.  This  code,  which  allows  2-d  or  3-d 
primary  disturbances  and  an  arbitrary  number  of  secondary  modes  of  different  types,  has 
been  verified  against  existing  data  for  incompressible  boundary  layers  and  against  direct 
numerical  simulations  at  M <*,  =  1.6  and  4.5. 

The  influence  of  the  (2-d)  primary  mode  amplitude  A  on  the  secondary  instability  char¬ 
acteristics  has  been  studied  at  four  Mach  numbers:  Moo  —  0, 0.8, 1.6, 4.5.  We  found  that  a 
small  but  finite  primary  disturbance  amplitude  can  trigger  a  broad-band  of  highly  unstable 
3-d  secondary  disturbances  whose  growth  rates  are  large  and  increase  with  increasing  A. 
Increasing  the  free  stream  Mach  number  not  only  pushes  the  disturbances  away  from  the 
wall,  but  also  confines  them  to  a  much  smaller  region,  centered  near  the  critical  point. 

The  influence  of  the  primary  disturbance  skewness  9  (with  respect  to  the  mean  flow 
direction)  on  secondary  instability  characteristics  at  =  1.6  is  first  examined  for  9  = 
0°,  5°,  10°,  15°.  The  unstable  subharmonic  and  the  fundamental  disturbances  propagate  with 
the  same  phase  velocity  as  the  2-d  primary,  whereas  their  phase  velocities  are  different  from 
the  phase  velocity  of  the  3-d  primary.  For  6  =  10°  and  8  =  15°,  there  are  two  maxima 
in  the  curve  of  the  growth  rate,  aT,  versus  the  wavenumber  in  the  direction  parallel  to  the 
phase  front  of  the  primary  wave  (i.e.,  /32).  The  strongest  subharmonic  modes  occur  when 
the  primary  disturbance  propagates  in  the  mean  flow  direction,  while  the  most  unstable 
fundamental  modes  occur  when  the  primary  disturbance  is  inclined  at  about  10°  to  the 
mean  flow  direction.  When  the  subharmonic  instability  is  examined  at  a  relatively  high 
value  of  8  (namely,  8  —  45°),  it  is  found  the  first  five  Fourier  modes  (i.e.  N  =  5)  are  required 
to  capture  the  secondary  growth  rate  correctly  to  4  significant  digits. 

At  Moo  —  4.5,  there  are  two  types  of  primary  instability  modes:  first  and  second  mode. 
At  this  Mach  number,  the  subharmonic  and  fundamental  instabilities  which  originate  from 
both  the  first  and  the  second  mode  waves  are  examined.  We  found  that  the  subharmonic 
instabilities  which  arise  from  the  second  mode  disturbances  are  the  dominant  instabilities. 
When  the  primary  wave  amplitude  increases  from  2  %  to  6  %,  the  growth  rate  of  the  most 
amplified  subharmonic  of  a  second  mode  disturbance  is  augmented  by  about  136%. 

In  summary,  we  believe  that  secondary  instability  is  a  viable  boundary-layer  transition 
mechanism  at  =  0.8, 1.6  and  4.5.  Further,  the  high  frequency,  acoustic  type  second 
mode  subharmonic  (secondary)  disturbances  will  prevail  at  Ma 0  =  4.5. 
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Fig. 5:  Comparison  of  theory  with  direct  numerical  simulation. 
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Growth  rate  as  a  function  of  primary  wave  amplitude.  The  results  are  for  M ^  — 
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Fig. 17:  Subharmonic  eigenfunction  at  M, =  0.0.  The  parameters  are  a  =  0.29,  6  =  0°, 
Tx  =  520°jR,  Re  =  1675,  and  A  =  1.5%. 
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Fig. 18:  Subharmonic  eigenfunction  at  =  1.6.  The  parameters  are  a  = 
(32  =  0.4,  T =  520 °R,  Re  =  1675,  and  A  =  1.5%. 
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Fig. 25:  Normalized  streamwise  velocity  components  of  the  subharmonic  instability  modes 
at  Moo  =  1.6,  a  =  0.24,  0  =  0°,  ft  =  0.4,  7^  =  520°#,  Re  =  1675,  and  A  =  1.5%. 
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